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Abstract 

The dynamical effects of ground state correlations for excitation energies 
and transition strengths near the superfluid phase transition are studied in 
the soluble two level pairing model, in the context of the particle-particle self 
consistent Random Phase Approximation (SCRPA). Exact results are well 
reproduced across the transition region, beyond the collapse of the standard 
particle-particle Random Phase Approximation. The effects of two-body cor- 
relation in the SCRPA are displayed explicitly. 



PACS numbers: 21.60.Jz, 21.60.Fw, 21.60 Ev 



1 



An extension of the self-consistent Random Phase Approximation (SCRPA) to the 
particle-particle case (pp-SCRPA) has been recently presented by Dukelsky, Ropke and 
Schuck |1|]. The main concern there is with differences between ground state correlations due 
to pairing vibrations and Bruekner-Hartree-Fock ground state correlations. In this report 
we use the soluble schematic model used in [|TJ to study different aspects of the pp-SCRPA. 
Namely, we examine the behavior of the system as we approach the transition point from 
the normal to the superfluid phase and investigate the importance of the two-body ground 
state correlations for excitation energies and transition strengths in this domain. 

The model is the standard two-level pairing model whose Hamiltonian is given by 

H=^J2aN a - g nJ2AtA al (1) 

cr o-.cr' 

where a = ±1 labels the upper and lower levels which have equal pair degeneracy Q = j + 1/2 
and N a , are respectively the number and pair operators in the level a, i.e. 

m V » <• " m 

For each level a the operators 

sf = Va4+, s& = s+v + , s^ = ^(N a -n) (3) 

satisfy an SU(2) algebra. 

The addition, P + , and removal, R + , phonon creation operators of the pp-SCRPA for the 
model are J| 

p+ = \A+ - nA± x , R + = pA-x - tAx (4) 

and the pp-SCRPA ground state | 0) is defined as the vacuum of these phonons, i.e. 

P | 0) = R | 0) = 0. (5) 

Addition and removal excitations are normalized as 

(0|[P,P+]|0> = 1, (0|[i?,i?+]|0) = l, (0|[P+,P+]|0> = (6) 
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which gives, using Eqs. (|J), 

■"-(1-§)|0> + ^<0|(1-%1)|0> = 1 



A 



^0|(1-§)|0> + A-<0|(1-^)|0> = -1 
tA(0| (l - ^) |0> + W (0| (l - ^f) |0) = 0. 

The pp-SCRPA equations for the present model can be written in the form 
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with the usual definition of the symmetrized double commutator 

[A,B,C} = ^([A, [B,C]] + [[A,B],C]). 

Equation([|) has two solutions. We associate the solution of positive norm (normalized 
to +1) with the addition phonon, A + = A, /i + = /i and the corresponding eigenfrequency 
with uj + = Eq(A + 2) — Eq(A). The negative norm solution (normalized to -1) we associate 
with the removal phonon A_ = r, = p and the corresponding eigenfrequency with 
oj_ = Eq(A) — E (A — 2) pj. This is the physical meaning of the first two Eqs. (0). The 
last one expresses the orthogonality of these two solutions. 

^From Eqs.(0) we get r and p in terms of A and /i and, choosing phases to agree with 
ref. [Ill, one has 



p+ = \A+ - fiAtt 
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The model Hamiltonian ([!]) commutes with , ^ and Sq + Sq , where = 
- (S+^S^ + S^S+^ + So^ 2 . In particular, the ground states of even nuclei are defined 
in the SU(2) ® SU(2) irreducible representation (—, — ^. A basis of states in this subspace 
is given by 



mi; m 2 ) 
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The pp-SCRPA ground state of the N = 2Q system can therefore be expanded as 



n 
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where m is the number of particle pairs coupled to J = in the upper-level. 

The condition that | 0) is the vacuum of the addition phonon, Eq. ([5]) gives the recursion 
relation 



C — p 

^m-^—m » ^m— 1; — (m— 1) 



which allows for writing the normalized ground state explicitly as [1 

2 n 
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Apparently we have a problem here, since the pp-SCRPA ground state has been deter- 
mined by just one of Eqs. @ and we still have the additional condition that | 0) is also 
supposed to be the vacuum of the removal phonon. However, the second Eq. (||) gives for 
the coefficients of | 0) the condition 
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This is consistent with Eq. ( |TTD only if 

(0 | ALi - O | 0) = (0 | O - Ni | 0) 



(14) 



and that this is true follows from 



(0 I Nt + ALi - 2ft | 0) = 



which can be interpreted as saying that the number of particles in the level +1 is equal 
to the number of holes in the level —1. Note that the above discussion uses the condition 
iV = 2ft (or Sq 1 -* + = 0) in an essential way. 
Defining the normalized amplitudes 



A = AW{0 | 1 



ft 



0), /J = /iW(0| 1 
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with A — jj, — 1, we can write the pp-SCRPA equations as 



A B 

B* C 
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where A, B and C are now given by Eqs. 
defined as 



in terms of the normalized operators A 1 , A_ l 
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Evaluating the double commutators one obtains 



A = e-g 
B = g 
C = e-g 
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+ g(0 | 2A_ 1 A_ 1 + AiA_! + AZ X A X | 0). 



(16) 
(17) 
(18) 



In order to calculate the expectation values appearing in these expressions we invert Eqs. 
(|j) to find A^ = XP + — ~pR and A^i = ~pP + — XR. /,From this we then obtain, using the 
vacuum condition <§), (0 | A~^Ai | 0) = 7J 2 , (0 | A^A^ | 0) = A 2 and (0 | AtxAx | 0) = A/Z 
|Ij]. The matrices A, B and C are thus functions of the pp-RPA amplitudes A and JZ, of the 
one-body densities (0 | iVi | 0), (0 | ALi | 0) and of the two-body densities (0 | N 2 — iVi | 0>, 
(0 | iV 2 ! - ALi | 0) and (0|iVi ALi|0). Their final expressions are 
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Once we have the explicit expression (^) for the correlated ground state | 0), we can 
calculate the one- and two-body densities and solve the non-linear pp-SCRPA equations, 



Eqs. Q15[) and (|19|)-(piD, for the pp-RPA amplitudes A and ~p. Results of such a calculation 
will be presented further on. Note that in the limit of standard pp-RPA, | 0) is replaced by 
the uncorrelated ground state and Eq. (p~5|) reduces to the ordinary equation 
/ 
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e — gfl —g& 
-gQ e + 2g — gVt 

This illustrates the fact that, in practice, solving the SCRPA equations is rather de- 
manding from a computational point of view, and in fact this has not been done so far in 
realistic cases. Several authors have therefore suggested alternate schemes to derive simpler 
equations which are more easily amenable to numerical treatment ||, [|J- ]/]]. In essence, 
the simplifying hypothesis consists in neglecting the two-body correlations in the pp-SCRPA 
equations of motion (|15|). These approximations are generically called renormalized RPA 
(RRPA). The pp-RRPA simplified version of the ingredients involved in Eqs. (|16|)-(|T8|) reads 
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and A, B and C become respectively 



^- 9 <^>- 9 «<o, (l -£),o> 

B~- g n(o \ (i-^) |o> 

^ £ + 29 - 9 i^i)- 9 , ( o| (,-*), o). 

The pp-RRPA equation is thus analogous to the pp-RPA equation fl22|) , except that 
now g is renormalized by ^1 — - — ■ ^ ■ — -J and the particle and hole energies are shifted 

respectively by — g- — ■ ^ ■ — - and g- — ■ ^ ■ — -. In realistic cases || we have to make 
further approximations in the calculation of the one-body densities to avoid having to solve 
the vacuum condition, Eq. (|5]). For the calculations reported below we calculate (0 | Ni | 0) 
using the expression Eq. (|T2| ) for the phonon vacuum | 0). The difference between our 
pp-SCRPA and pp-RRPA results are thus solely due to neglecting two-body correlations in 
the pp-SCRPA equations of motion. 

We next present numerical results obtained for the two-level pairing model. They cor- 
respond to a case in which twenty particles distribute themselves in two Q = 10 levels 
separated by e = 2. In Fig. 1 we plot the energies uo + and uj- as functions of the coupling 
strength g. The pp-RPA collapses at g cr n = — — -, at which point the corresponding ener- 
gies become equal, u + = cu_ = — g crit . Since u + and u_ have the physical interpretation of 
(minus) two-nucleon removal energies respectively from the closed-shell-plus-two and from 
the closed-shell nuclei, this equality signals the breakdown of shell structure. On the other 
hand, neither the SCRPA nor the RRPA collapse, but instead follow closely the exact values 
in the neighborhood of the phase transition point. As g increases, however, the agreement 
with the exact energies deteriorates. In particular, the pp-SCRPA, and even more markedly 
the pp-RRPA energies approach values that differ from the exact strong coupling limits 
uj + — > and — > —2g for large values of g. In Fig. 2 we plot the occupation of the upper 
level in the ground state of the closed shell nuclei. Again we see that in the neighborhood of 
the phase transition point both approximations yield values close to the exact ones, whereas 
the pp-RPA value diverges. However, when g increases, the pp-SCRPA value lags consid- 
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erably behind the exact one, which agrees well with the pp-RRPA result almost up to the 
strong coupling limit of (0 | Ni | 0} — > Q. Fig. 3 shows, also as a function of g, the square 
of the two-particle transfer matrix element to the upper level: 

\ (N + 2\ A~f \ N) p = Wl-<°'^ |0 >). 

We find again the same behavior, with agreement near the phase transition point and dis- 

n + i 

agreement when g increases approaching the strong coupling limit — - — . 

In conclusion, we find that correlations introduced by the pp-SCRPA or by the pp- 
RRPA allow for following the system across of the phase transition point associated with 
the collapse of the ordinary RPA. However, it should be stressed that the pp-RRPA, which 
is in fact an approximate treatment of the full pp-SCRPA performs even better than the 
full theory for some observables. The pp-RRPA values for the occupancy of the upper level 
agree with the exact values all the way up to the strong coupling regime, and pp-RRPA 
values for the two-particle transfer matrix elements also perform somewhat better than the 
pp-SCPRA values. Only for the excitation energies are the pp-SCRPA results better than 
the corresponding pp-RRPA results. A possible perspective concerning these results could 
be formed by recalling the well known fact that the standard RPA tends to overestimate 
collectivity. In the present schematic model, on the other hand, the full pp-SCRPA seems 
to underestimate it, whereas the ommision of two-body correlations in the pp-RRPA gives 
improved occupancy and transfer results. In general, when g increases towards the strong 
coupling limit, pp-SCRPA and pp-RRPA correlations are not adequate to reproduce exact 
values, indicating that we should consider the system as a pairing-deformed (superfluid) 
system || . This is consistent with the fact that the pair-coupling model || is asymptotically 
exact in the strong coupling limit. 
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Figure captions 

Figure 1. Excitation energies u + (upper points and curves) and cj_ as functions of the 
coupling strength g for the two-level pairing model with e = 2, Q = 10 and N = 20. 

Figure 2. Number of particles in the upper level for the ground state of the two-level 
pairing model. Same parameters as for Fig. 1. 

Figure 3. Absolute square of the two particle addition matrix element to the ground 
state of the two-level pairing model. Same parameters as for Fig. 1. 
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de Passos et al.: Self Consistent... Fig. 1 
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de Passos et al.: Self Consistent... Fig. 2 




de Passos et al.: Self Consistent... Fig. 3 
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